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<N ■ ABSTRACT 
U 

' We report the angular diameters of two Mira variables (U Ari and Z Sco), three 

semi-regular and irregular variables (SW Vir, r\ Gem and [i Gem) and a supergiant 
semi-regular variable (TV Gem) by lunar occultation observations in the near-infrared 
broad K-band (2.2 \xm). Lunar occultations of n Gem and Gem were also observed 
for the first time simultaneously in both K and L' bands yielding angular diameters 
£NJ , at 2.2 and 3.8 /im. Effective temperatures and linear radii are also derived for all the 

observed sources and compared with earlier measurements. The mode of pulsation of 
CD ' both Mira and SR sources in our sample is discussed. 
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C3 . 1. Introduction 



Asymptotic giant branch (AGB) stars are in the last stage of stellar evolution before turning 
into planetary nebulae and are generally surrounded by circumstellar matter due to their large mass- 



loss rates (~ 10~ 6 M yr _1 ). AGB stars include classical Mira variables (visual amplitudes > 2.5 
mag ; periods, 100 - 1000 days), semi-regular variables SRa (visual amplitude < 2.5 mag ; periods 
35 - 1200 days), semi-regular variables SRb (amplitude < 2.5 mag ; with poorly defined periods), 
and irregular variables Lb type (amplitude is small, no definite periods), as well as supergiant 
semi-regular variables SRc type (Smith et al. 2002). The evolutionary connection between Mira 
and SR group is not clear, but SR variables are often considered to be Mira progenitors (Bedding 
& Zijlstra 1998). 

Multiwavelength measurements of angular sizes of Miras and Semi-Regular variables at differ- 
ent phases of their pulsation cycle provide a direct means to understand their atmospheric extension 
and pulsation properties. The high mass-loss and relatively low surface temperature of the evolved 
stars provide a habitable zone for several molecules like TiO, VO, H2O & CO etc. in their extended 
atmospheres. Large opacities of these atmospheric molecules in some particular bands mask the 



- 2 - 



dominant photospheric continuum radiation. The emergent observed radiation is thus contami- 
nated by radiation from the relatively cold atmospheric layers. Consequently photospheric size 
measurements are effected in the different filter bands which has been known for sometime (e.g. 
Labeyrie et al. 1977; Quirrenbach et al. 1993; Haniff, Scholz & Tuthill 1995). 

Recently there have been many high quality interferometric measurements of angular diameters 
of Miras at near infrared wavelengths (e.g. Mennesson et al. 2002; Woodruff et al. 2004; Perrin et 
al. 2004; Millan-Gabet et al. 2005; Fedele et al. 2005). Mennesson et al. (2002) found that L'- 
band diameter of several oxygen rich Miras were much larger (25% to 100%) than those measured in 
the broad K-band and proposed as an explanation a simple empirical model of a central stellar disk 
surrounded by a optically thin gaseous shell. Observed variation of angular size with wavelength can 
be interpreted in terms of transparency of the optically thin shell varying with wavelength. Perrin 
et al. (2004) have observed several Miras in narrow bands around 2.2 /xm and find systematically 
larger diameters in bands contaminated by water vapor or CO. Millan-Gabet et al. (2005) report 
a systematic increase of angular size with wavelength (~ 25%) from J to H to K' from a study of 
23 Miras involving simultaneous measurements in JHK'. The 11 /jm interferometric observations 
of three Mira variables (Weiner, Hale & Townes 2003a; Weine, Hale & Townes 2003b) showed 
the diameters at 11 fim larger by a factor of ~2 than those measured in the K-band. The increase 
of apparent diameter from the near-infrared toward longer wavelengths seems to be a common 
phenomenon in Mira and late M-type semiregular variables. The dispersion in angular sizes from 
near-infrared to mid-infrared band are well represented by the modeling of interferometric data by 
inclusion of water shell surrounding the Mira variables (Weiner 2004; Ohnaka 2004; Schuller et 
al. 2004). ISO and ground-based spectroscopy observations also point to the warm molecular shell 
of H2O surrounding the Mira variables (Hinkle & Barnes 1979; Tsuji et al. 1997; Yamamura, de 
Jong & Cami 1999; Matsurra et al. 2002; Tej et al. 2003b). Several theoretical models have been 
developed to understand the dynamic atmosphere of Mira and non-Mira M stars (Bessell, Scholz 
& Wood 1996; Hofmann & Scholz 1998; Hofmann, Scholz & Wood 1998; Hofner et al. 1998, 
2003; Woitke et al. 1999). Recent results of interferometric observations have been predicted well 
by considering molecular contamination effects on theoretical models (Jacob &: Scholz 2002; Tej, 
Lancon & Scholz 2003a). 

The question of pulsation modes for Mira stars is a complicated one. Theoretical considerations 
generally suggest a fundamental mode as it is difficult to reproduce in the first overtone, the large 
velocity amplitude encountered in Miras (Bessell, Scholz &: Wood 1996). The quantity commonly 
used in stellar pulsation modelling is the Rosseland radius of the hypothetical parent star of the 
Mira variable which does not pulsate.lt is not an observed quantity but related to the intensity 
distribution on the stellar disk (Scholz 2003, and references therein). It is also a phase and possibly 
cycle dependent value. Jacob & Scholz (2002) showed that transforming a measured diameter into 
a Rosseland value may be difficult or impossible due to molecular contamination in even standard 
near continuum bandpasses. The variability of diameter with Mira phase and cycle as well as 
wavelength further complicates the issue. Angular diameter that are reported from observed data 
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are obtained by fitting visibilities (from interferometric data) or occultation lightcurves (from lunar 
occultation) to a well denned artificial center to limb variation (CLV) like a uniform disk (UD), 
a fully darkened disk (FDD) or a Gaussian intensity distribution. It is difficult to obtain the 
actual CLV from the observed data. Another difficulty in determining pulsation modes is the large 
uncertainty in distance measurements required to convert angular to linear diameters. While earlier 
work (Haniff, Scholz & Tuthill 1995; van Leeuwen et al. 1997; Tej et al. 1999; van Belle, Thompson 
& Creech-Eakman 2002) suggested first overtone pulsations in some Miras, recent results based 
on interferometric diameter measurements (Perrin et al. 2004; Woodruff et al. 2004; Fedele et al. 
2005) which take into into account molecular contamination in the band passes (Jacob & Scholz 
2002; Mennesson et al. 2002; Tej, Lancon & Scholz 2003a; Ireland et al. 2004a,b,c; Ohnaka 2004) 
point to a fundamental mode of pulsation in Miras. 

In this paper we present new angular size measurements of six evolved variables by lunar occul- 
tation (LO) in the near-infrared. These include two Miras (U Ari and Z Sco) and four semiregulars 
(SW Vir, i] Gem, /x Gem and TV Gem) including the supergiant TV Gem. Our occultation re- 
sults from another Mira U Ori are also included for comparison although it's angular diameter has 
been reported earlier (Mondal & Chandrasekhar 2004). We report the first simultaneous angular 
diameter measurements at 2.2 and 3.8 /im for two SR variables (rj Gem and \x Gem), rj Gem 
underwent two LO events and in both cases simultaneous K and L' angular diameters could be de- 
termined. Uniform disk angular diameters are derived from our occultation data. Bolometric fluxes 
are estimated from photometry and effective temperature are calculated. Distances estimated to 
the sources are considered and linear radii are derived. The position of these evolved variables is 
plotted in the period-linear radius diagram and their mode of pulsation discussed. 



2. Observations and data analysis 

Three sources (U Ari, rj Gem and \i Gem) were observed simultaneously in broad K (2.2/0.40 
\xm) and L' (3.8/0.60 (im) filter-bands while Z Sco, SW Vir, TV Gem and U Ori occultations were 
recorded in the K-band only. All observations were made using two channel high speed photometer 
installed on the 1.2m telescope at Mount Abu, India. The details of the instrument can be found 
elsewhere (Mondal et al. 1999; Mondal, Chandrasekhar & Kikani 2002). The sampling time of 
the light curves was 2 milliseconds (ms) except SW Vir which was sampled at 1 ms. rj Gem was 
observed twice during an interval of 2 months. All events expect Z Sco were recorded under good 
sky conditions, Z Sco event was recorded through thin clouds. The details of occultation events are 
listed in Table 1. 

The observed light curves containing modified Fresnel diffraction fringes are modeled to get the 
uniform disk angular diameter of the stellar sources. The model fitting of the lunar occultation light 
curve involves the x 2 -iffimmization technique to obtain the best estimation of the five parameters 
: (i) the geometric time of occultation; (ii) the stellar signal; (iii) the sky-background; (iv) the 
velocity component of the moon in the direction of occultation and (v) the uniform disk angular 
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diameter. The analysis procedures are based on the standard non-linear least- square (NLS) method 
introduced by Nather & McCants (1970). The point-source Fresnel diffraction pattern modulated 
by the finite spectral bandwidth of the system, the finite telescope aperture, the instrument time 
response and the extended angular size of the source are taken into consideration for fitting of the 
above-mentioned parameters. The resolution limit of the LO technique experimentally determined 
by studying occultations of a number of bright point sources, is ~ 2 mas (Chandrasekhar 1999). 

3. Results and Discussions 

Individual source parameters are listed in Table 2. Mass-loss rates and outflowing velocities 
derived generally from CO and SiO line measurements are also listed in Table 2. The difficulty of 
obtaining good distance estimate is discussed later in section 3.3 and adopted distances to sources 
are given in Table 2. 

3.1. Angular diameter measurements 

3.1.1. U Ari 

U Ari is an oxygen-rich Mira variable of period 371 days with a spectral type M4-9.5 IIIc 
(Keenan & McNeil 1989). 

Lunar occultation of U Ari was observed simultaneously in the K and L' bands close to min- 
imum phase ( phase 0.57). A good occultation trace has been recorded in K-band (Fig. 1). The 
signal to noise ratio (S/N) limited by atmospheric scintillation noise in K-band is about 50. The 
L'-band profile is noisy and has not been considered for analysis. 

We obtain a UD angular diameter of 7.3±0.3 mas in the K-band. The UD model fit to the 
K-band light curve is shown in Figure 1. An earlier LO measurement in the H-band gave a UD 
diameter of 6.11±0.34 mas at phase 0.49 (Ridgway et al. 1979). There is a significant difference 
of about 20% in the two UD values which have been measured at nearly the same phase near the 
minimum. From recent simultaneous size measurements of 23 Miras at JHK'-bands, the larger size 
in K-band compared to H-band appears to be a common phenomenon in Mira stars (Millan-Gabet 
et al. 2005). The apparent size variations from the near to mid- infrared wavelengths are reasonably 
well modeled by inclusion of warm (1500-2000 K) H2O shell within a few stellar radii (Mennesson et 
al. 2002; Jacob & Scholz 2002; Tej, Lancon & Scholz 2003a; Weiner 2004; Ohnaka 2004; Perrin 
et al. 2004). Furthermore, it is to be noted that as both observing (K and H-band) phases arc 
close to the minimum, the molecular contaminations are expected to be more prominent compared 
to the maximum phase and also evident from spectroscopy (Tej et al. 2003b). Hence our simple fit 
of UD diameter may not represent the true continuum size near the minimum of the Mira phase. 
In case of interferometric observations, the observed visibilities are compared successfully to the 
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theoretical models of visibilities with inclusion of a thin surrounding H2O shell (Perrin et al. 2004). 
However, in case of LO such modeling efforts have not yet been carried out. 

3.1.2. Z Sco 

Z Sco is an oxygen-rich Mira variable of period 352 days with a spectral type M4/5 Hie. The 
maser lines of CO and SiO are not detected in the circumstellar atmosphere of Z Sco (Young 1995; 
Cho, Kaifu & Ukita 1996) and hence there is no mass-loss estimate and outflow velocity of the 
source. No dust signature is present in IRAS LRS spectra (Sloan & Price 1998). 

The occultation of Z Sco was recorded in the K band at phase 0.26. The sky condition during 
observations was poor (thin passing clouds) but it was possible to record the event. Fringe distortion 
is evident in the occultation trace (Figure 2) but nevertheless four fringes are recorded. The light 
curve is fitted with the UD model and a varying background using 5th order Legendre polynomial 
shown in Figure 2. 

We first derive the UD angular diameter of 3.8±1.0 mas. The error is large due to poor quality 
of the light curve. The diameter of Z Sco has never been measured before by a direct technique. 
Using the empirical relation of van Belle (1999), based on (V-K) magnitude, we derive an angular 
diameter of 5.1±1.3 mas using the visual magnitude of 11.6 from the AAVSO database (Mathei 
2004). We measure the K magnitude of the source to be 1.33±0.1 at the time of our occultation 
observations. 

3.1.3. U Ori 

U Ori is an oxygen-rich Mira variable of pulsation period 371 days with a spectral type 
M6-M9.5 Hie (Keenan k McNeil 1989). 

The occultation of U Ori was recorded in the K-band at phase 0.33. We obtain the UD 
diameter of 11.90±0.30 mas. We had earlier studied this source from the point of view of the 
spatial asymmetry. From the comparison of near simultaneous lunar occultation observations on 
U Ori at the same wavelength (K-band) and at different position angles (PA) (75° and 136°), from 
two observations we found the evidence of asymmetry in its atmosphere; the source appears to be 
elongated at PA of about 70° (Mondal & Chandrasekhar 2004) . The asymmetric extension is also 
found in OH maser observations on U Ori (Chapman, Cohen & Saikia 1991). 



3.1.4. SW Vir 

SW Vir is an oxygen-rich semi-regular pulsating variable of SRb type having a pulsation 
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period of 150 days and spectral type of M7 III (Lebzelter & Hron 1999). Kiss et al. (1999) 
identified it as a triply periodic variable (1700, 164 and 154 days). Here we have adopted the 
period of 150 days that is widely accepted in the literature. 

The lunar occultation of SW Vir (M7 III) was recorded in the K-band only. The S/N of data 
is limited by the atmospheric scintillation. We derive the UD diameter of 15.9±0.6 mas (Table 3). 
The best-fit UD model with observed data-points are shown in Figure 5. 

Three angular diameter measurements of SW Vir was reported previously in the near-IR 
region by LO and the reported values are in the range of 16.11 to 16.82 mas (Ridgway et al. 1982; 
Schmidtke et al. 1986). Details of those measurements are listed in the Table 4. Recently using 
LBI the reported values of angular diameters are 16.24±0.06 mas in the K' band and 22.88±0.33 
mas in the L/ band (Mennesson et al. 2002). The molecular contamination in SW Vir is prominent 
from K and L' size measurements. Our measurement in the broad K-band is consistent with 
previous measurements. From previously reported measurements (in Table 4) we also note that 
phase-dependent variations are not observed in SW Vir. 



3.1.5. i] Gem 

7] Gem is an oxygen-rich semi-regular variable (SRa) having a spectral type of M2.5 III 
(Keenan & McNeil 1989). The variability is classified as SRa with small visual amplitude of 0.75 
mag. Two periods, 233 and 20 days, are found from photometric observations (Percy & Wilson 
2001). Most of the earlier observations report the longer period, 233 days. 

Two lunar occultations of r] Gem (M2.5 III) have been observed on 08 January 2001 and 4 
March 2001 simultaneously in the K and L' bands. The best-fit UD model curves with observed 
data points are shown in Figure 3. We derive the UD angular diameters of 12.7±0.3 mas at K and 
12.7±1.0 mas at L' for 08 Jan 2001 and 12.8±0.3 mas at K and 12.8±2.0 mas at L' for 04 Mar 2001 
(Table 3). 

We thus have two good sets of observations in both K and L' bands separated by a two months 
interval. In this period a good K-band measurement of angular diameter has been made (Richichi 
& Calamai 2003) which is in good agreement with our value. From our two-epoch, two-wavelength 
observations we can conclude that there is no detectable variation of angular size with phase in 
case of r\ Gem. 

Previous angular diameter measurements of 77 Gem are listed in Table 4. The source has a 
well determined angular diameter from the wavelength range 0.55 to 2.2 /j,m. The optical (0.55 //m 
and 0.80 /urn) UD diameters are 11.43±0.55 and 10.91±0.11 mas respectively (Mozurkewich et al. 
2003). At 0.712 jiva. (in strong TiO band) and at 0.754 fim (in adjacent continuum), UD diameters 
are 11.75=L0.27 mas and 10.70=1=0. 15 mas respectively (Quirrenbach et al. 1993). These optical 
diameters are not different from continuum diameters and are slightly lower than our measured IR 
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diameters. 

r\ Gem is also identified as a spectroscopic binary. The spectral type of the companion is 
identified as GO III with visual magnitude of 11.3 and the separation of 0.9 to 1.08 arcsec from 
primary at position angle (PA) of 29° (Phillips et al. 1980; Baize 1980). The Hipparcos catalog 
shows a binary separation of 1.7 arcsec at PA of 261° (Perryman et al. 1997). The variation of 
separation in several observations is attributed to the ellipticity of the orbit. No binarity signature 
is detected in any of our LO light curves and it was also undetected from previous LO observations. 
The brightness ratio between the primary and secondary component is estimated ~ 1:1600 in the 
visual band. The K magnitude of the companion would be ~ 13 mag which is well below the limit 
of our detection. 



3.1.6. n Gem 

fj, Gem is an oxygen-rich semi-regular variable of Lb type. The source has a spectral type of 
M3 III (Keenan k McNeil 1989), a period of 27 days (Percy & Wilson 2001). 

The lunar occultation of fi Gem (M3 III) was also observed simultaneously in the K and L'- 
bands. We derive the UD angular diameters of 13.7±0.5 and 14.8±1.0 mas in the K and L'-band 
respectively (Table 3). The model-fit light curves along with observed data points are shown in 
Figure 4. 

There are thirteen observations of lunar occultation in the wavelength range 0.4 to 0.82 /im 
listed in the catalog of White & Feierman (1987) and the mean UD value in that wavelength 
range is 13.06±0.42 mas. Recently UD angular sizes in the optical bands from interferometric 
observations are 13.98±0.14 mas (at 0.80 fim) and 13.48±0.19 mas (at 0.55 /um) (Mozurkewich et 
al. 2003). UD sizes at TiO absorption band ( 0.712 /xm) and nearby continuum (0.754 /im) are 
13.97±0.28 mas and 13.50±0.13 mas respectively (Quirrenbach et al. 1993). UD value at K-band 
is 13.50±0.15 mas (Di Benedetto & Rabbia 1987). Mira-like enlargement (a factor of ~2) at TiO 
band compared to the adjacent continuum has not been noted in fj, Gem. Some of the previous 
measurements are listed in the Table 4. 

Considering all available measurements including our own it appears that the UD diameter of 
H Gem has not shown any substantial variation from optical to near-IR over many years. Our UD 
angular diameters also show no significant variation between K and L' bands within the errors of 
measurements. 

For non-Mira stars up to at least spectral type M4, it has been shown that observed and 
model-predicted visibility curves do not differ significantly from the UD profile (Wittkowski et al. 
2001, 2004). This is also realized from theoretical models of non-Mira M giants (Hofmann & Scholz 
1998) as against the Mira model (Hofmann, Scholz & Wood 1998). In case of non-Mira stars (77 
Gem, /j, Gem and TV Gem but not SW Vir) in our sample, UD size in the K-band may be close 
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the true continuum diameter consistent with previous observations (Table 4). Although warm H2O 
is detected in the spectra of K and early M giants (Tsuji 2001), for rj Gem and n Gem we find 
that there is no change between K and L' diameters unlike in Miras or SW Vir (M7 III). Molecular 
layers if at all present in the atmospheres of i] Gem and fi Gem must have column densities too 
low to affect L band diameters unlike in Miras. 



3.1.7. TV Gem 

TV Gem is an oxygen-rich supergiant semi-regular variable (SRc). The spectral type is Ml-0 
lab (Keenan & McNeil 1989). The visual magnitude varies from 7.0 - 7.8 over the pulsation period 
of 182 days (Kukarkin et al. 1969). The distance we have adopted here, is 1200± 300 pc based on 
interstellar extinction towards Gem OBI association (Underhill 1984). 

The lunar occultation of supergiant TV Gem was recorded on 14 Nov 2000 in the K-band under 
clear sky conditions. The S/N is high (~100) and is limited by atmospheric scintillations. The 
UD model fit to the light curve with usual five free parameters (in section 2) was not completely 
satisfactory. A better fit to the data is obtained by including a star plus shell model rather 
than a single star model which is shown in the residuals of fits in lower panel of Figure 6. The 
JHK photometric observations on 16 Nov of 2000 yielded the magnitudes 2.31±0.05, 1.39±0.05 & 
1.16±0.06 in J, H and K respectively. 

Angular diameter measurements of TV Gem by lunar occultation have been reported several 
times and are listed in Table 4. Earlier lunar occultations reported UD values of 5.31±0.91 mas 
in the optical region (Radick, Henry & Sherlin 1984), 4.9±0.3 mas in the K band (Ragland, 
Chandrasekhar & Ashok 1997) and 4.46±0.07mas again in the K band (Richichi et al. 1998). 
From LO observations in the K-band Ragland, Chandrasekhar & Ashok (1997) had reported the 
double shell structure of TV Gem like another supergiant a Ori (Danchi et al. 1994). The inner 
dust shell was estimated to be at 20±5 R*. The outer shell was estimated to be at ~ 500R* based 
on LRS spectra IRAS and IRAS photometry (12, 25, 60 jiva). They found the shell contribution 
in the K-band to be ~3%. We measure the UD angular diameter of the source 4.8±0.2 mas. We 
estimate the dust shell to be at 13±5 R*. The shell contribution to the K-band flux is ~5% . These 
results are consistent with earlier measurements by Ragland, Chandrasekhar & Ashok (1997). We 
estimate the effective temperature to be 3750±120 K again in good agreement with earlier values 
of Richichi et al. (1998). 

3.2. Bolometric flux and effective temperatures 

The bolometric fluxes are estimated by fitting a blackbody curve to available broad-band IR 
photometry measurements (JHKLM) compiled in the Infrared catalogue of Gezari, Pitts & Schmitz 
(1999) and 12, 25 Sz 60 fim IRAS PSC measurements including our JHK measurements in some 
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cases. In some cases (U Ari, U Ori and Z Sco) a two-temperature blackbody is required to best fit all 
observed points (1.25 -60 /im). Specifically an additional blackbody curve with cooler temperature 
(~500 K) fits the excess in IRAS flux. Such fits are shown in Figure 7. For rj Gem a single 
temperature blackbody curve is adequate for all observed fluxes (Fig. 7). The observed broad- 
band photometric magnitudes (JHKLL/) are converted to flux densities using the zero magnitude 
to flux density calibration established by Bessell, Castelli & Plez (1998). The blackbody flux is 
normalized with the observed flux in the K-band. In case of two temperatures blackbody fit to 
the photometric data, the star flux is fitted with relatively hotter temperature blackbody curve 
that is normalized with the observed K-band flux. To fit the excess in the observed infrared fluxes 
with additional relatively cooler temperature blackbody curve is further normalized with the star 
continuum subtracted flux at 25 /jm using the IRAS flux at that wavelength. By numerically 
integrating the single (or resultant of two temperatures) blackbody curve in the wavelength range 
0.4 to 100 |im, the bolometric flux is calculated. No reddening corrections were applied to estimate 
bolometric fluxes. These were deemed unnecessary, since typical magnitude of the corrections of 
our sample will be less than 0.05 mag in K-band. For example, the largest visual extinction was 
found in U Ori, A„ ~ 0.25 mag (Whitelock, Marang h Feast 2000a) and correspondingly ~ 
0.03 mag using the wavelength-dependent extinction relation, A& = O.llAt, established by Bessell, 
Castelli & Plez (1998). Our own infrared JHK photometric measurements of U Ari are used to 
estimate the bolometric flux at that particular phase (near minimum) while others are taken from 
Catchpole et al. (1979) at a similar phase. In case of SW Vir and TV Gem the bolometric fluxes 
were taken from the literature (Perrin et al. 1998; Ragland, Chandrasekhar & Ashok 1997). For 
other sources (U Ari, U Ori, Z Sco, ri Gem and r\ Gem ) we use our estimated bolometric fluxes 
from blackbody fits to calculate the effective temperature using the relation given below, 



where the bolometric flux F^ is in units of 10 -8 erg cm 2 sec -1 , the UD angular diameter <\> 
is in milliarcsec and the effective temperature T e jj in K. The typical error in bolometric flux is 
estimated to be about 15%. The effective temperature of the sources using our derived K-band 
UD diameters are listed in Table 5. For Mir as which are sources with extended atmospheres , T e ff 
refers to a specific layer. The characteristic reference level seems to be near optical depth unity in 
the near-IR continuum (Perrin et al. 2004). The effective temperature 2280±80 K of U Ari at 
phase 0.57 is substantially lower than the expected value of Mira stars (~3000 K) while the value 
2905±80 K of U Ori at phase 0.33 and 3120±420 K for Z Sco at phase 0.26 are consistent. 

We estimated the effective temperatures of 3450±125 K for rj Gem and 3675±140 K for /i 
Gem consistent with their spectral type. For SW Vir we estimate the effective temperature to be 
3060±130 K consistent with earlier measurements. 




(1) 
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3.3. 



Linear radii and mode of pulsation 



From interferometric observations and theoretical models it appears that the true continuum 
size estimate is difficult or impossible for Mira variables in presence of molecular contamination 
effects and phase cycle effects. Fedele et al. (2005) draw similar conclusions for R Leo. Other 
than wavelength and phase effects, determination of the mode for pulsation in Mira stars from 
period-radius relation is also constrained by large uncertainty in the distance to most of the sources 
as noted earlier in section 3 (Whitelock & Feast 2000b; van Leeuwen et al. 1997). Theoretical 
models (Bessell, Scholz &: Wood 1996; Hofmann, Scholz & Wood 1998) concluded on the funda- 
mental mode for these stars because it is difficult to produce the observed large velocity amplitude 
(Hinkle, Lebzelter & Scharlach 1997; Scholz & Wood 2000) from the first-overtone model. The 
comparison of theoretical pulsation models with MACHO observations of long-period variables in 
the LMC (Wood et al. 1999), pulsation velocities derived from Doppler line profiles (Scholz &; 
Wood 2000) also strongly indicate that Mira stars are fundamental mode pulsators. Earlier ob- 
served UD diameters for Mira stars at optical and near-infrared wavelengths without considering 
molecular contamination pointed to the first-overtone mode (Haniff, Scholz & Tuthill 1995; van 
Belle, Thompson & Creech-Eakman 2002). More recently considering molecular contamination on 
observed visibilities, Perrin et al. (2004) arrived at a lower size compared to earlier UD-fit size of 
six Miras and concluded that these are pulsating in fundamental mode. Woodruff et al. (2004) 
found that the interferometric visibilities curve of o Cet differs from UD profile and the data is 
better fitted with the fundamental model. 

SRs are separated from Mira stars by the shorter period and smaller amplitude, and they 
often show evidence of multiple periods (Bedding et al. 1998; Kiss et al. 1999). The mode of 
pulsation of Galactic Miras are well studied while little attention has been given to the pulsation 
studies on SRs. Feast (1996) found that SRs in globular clusters are pulsating in first overtone. 
From MACHO observations for LMC red variables, Wood et al. (1999) concluded that SRs can 
be pulsating in the 1st, 2nd or 3rd overtone, or the fundamental mode. Comparing observational 
and theoretical Q-values of 13 Galactic SRs, Percy & Parkes (1998) found that the majority are 
pulsating in the 1st or 2nd overtone while some are pulsating in the fundamental mode. 

Wood (1990) had suggested that in principle one combined equation for the position of low 
mass AGB stars (both Miras and semi-regular variables) can be used for comparing the observa- 
tional results. The standard pulsation equation is written as, 



where Q is a constant quantity (unit of days) that has distinct value for each mode of pulsation. 
Here R is actually the Rosseland radius of the non-pulsating star. Theoretical model predicted Q- 
values vary with period, mass and luminosity (Fox & Wood 1982). Typically the Q-value is pa 
0.105 (fundamental) and sa 0.04 (1st overtone) estimated from the theoretical models by Fox & 




(2) 
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Wood (1982). The Q- values used by Percy & Parkes (1998) for study of galactic SR variables are 
0.04 (1st overtone), 0.022 (2nd), 0.017(3rd) and 0.012 (4th) days from the models of Xiong, Deng 
& Cheng (1998). We have estimated Q-values for the samples considering mass of 1 M (Table 
6). Furthermore, following Ostile & Cox (1986) we have considered the following expressions for 
fundamental and overtone modes respectively, 



where P is the period in days. 

Masses of Miras of moderate (<400 days) period are reasonably well constrained ~ 1 M . 
Wyatt & Cahn (1983) estimated the main sequence masses of 124 Miras considering available data 
on radial velocity measurements. The Main-Sequence (MS) mass of progenitor of one Mira (U Ari) 
in our sample has been determined by Wyatt & Cahn (1983) to be 1.3M . Jura & Kleinmann 
(1992) suggested that the MS masses of Mira progenitors are in the range 0.8 - 2.0 M for Miras 
having the period less than 400 days. Theoretical models of Mira stars for masses of 1 and 1.2 Mq 
have predicted that geometric pulsation of continuum-forming layers is little affected by the mass 
difference (Ireland et al. 2004b). As semi-regular variables are progenitors of Miras, such mass 
considerations may be also applicable for them. In this mode analysis we have considered the mass 
range of 1.0 - 2.0 M . 

The linear radius is obtained from our K-band UD angular diameters and the adopted distance 
(Table 6) and superimposed on the theoretical model curves (Fig. 8). The errors shown on linear 
radii are mainly due to errors on distance estimates. 

There has been a great deal of discussion on Hipparcos parallax measurements of Miras in the 
literature(van Leeuwen et al. 1997; Whitelock &: Feast 2000b; Knapp et al. 2003). The revised 
Hipparcos parallaxes for U Ori and Z Sco have large errors (Table 6) and so we have not adopted 
them for distance estimation. Hipparcos parallaxes are excellent database for n Gem, r\ Gem and 
SW Vir, and we have adopted those parallaxes (Table 6). The adopted distances to all Miras in 
our sample have been estimated using the Period-Luminosity (PL) relationship given in Whitelock 
& Feast (2000b) (M K = - 3.47 logP + 0.84, M K is the absolute K magnitude and P the period) 
and the reddening corrected mean apparent K magnitude from Whitelock & Feast (2000b) (U 
Ori and Z Sco) and Feast (1996) (U Ari). The PL relations for Miras have been developed by 
several other authors (Feast 1996; van Leeuwen et al. 1997; Knapp et al. 2003) and those are 
consistent with the adopted relation. Considering uncertainties in apparent magnitudes (because 
of variability amplitude) and the dispersion in PL relations, we have considered errors in the Mira 
distances to be about 20%. 



logP = 1.86\og(R/R e ) - O.731og(M/M ) - 1.92 



(3) 



logP = 1.59\og(R/R Q ) - O.511og(M/M ) - 1.60 



(4) 



From Figure 8 and the Q- value in Table 6, it appears that the Mira star U Ari is pulsating in 
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the 1st overtone mode from both K and H-band measurements. However as noted in section 3.1.1 
observations in both K and H bands were near the minimum phase when molecular contamination 
effects are at a peak. Converting the measured LO UD diameter directly to a linear diameter would 
tend to overestimate the linear size. Till further measurements at other phases are available the 
result of overtone pulsation in U Ari must be treated with caution. 

The Mira star Z Sco is probably a fundamental mode pulsator but precise angular diameter 
determination was not possible because of noisy data. The mode of pulsation for U Ori appears to be 
a borderline case. We note that similar conclusions were drawn from earlier K-band interferometric 
observations of van Belle, Thompson & Creech-Eakman (2002). The minimum diameter for U Ori 
from recent interferometric observations in several narrow bands inside broad K-band (Perrin et 
al. 2004) however favours the fundamental mode. 

It is difficult to draw any conclusion on the SR variables SW Vir and n Gem from their 
position in Fig. 8. However, comparing observational (Table 6) and theoretical Q-values of Percy 
&; Wilson (2001), both SW Vir \x Gem may be the overtone candidates. The SRa rj Gem could be 
a fundamental mode pulsator of low mass but it has a second period of 20 days which complicates 
the issue. 



4. Summary and conclusions 

Our UD angular diameter of U Ari in the K-band shows a substantially larger value (~ 20%) 
compared to H-band observed earlier at nearly the same variability phase. Such an enhancement 
is consistent with a hot extended molecular layer close to the photosphere as suggested by Perrin 
et al. (2004). 

The supergiant TV Gem in the K band yields the angular size of 4.80±0.20 mas consistent 
with previous measurements. The dust shell around TV Gem is re-confirmed. We measure the 
dust shell size to be 13±5 R*. The effective temperature derived is 3750±120 K, consistent with 
the earlier value. 

We estimate linear radii of three Mira variables (U Ari, U Ori and Z Sco) from the K-band lunar 
occultation uniform disk angular diameters and distances derived from the PL relation. Comparing 
theoretical Period-Radius plots (Fig. 8) we find that U Ari is a first overtone pulsator, Z Sco is 
probably a fundamental mode pulsator while U Ori is a borderline case between fundamental and 
first overtone modes. A stronger conclusion regarding the pulsation mode of Miras from occultation 
observations could probably be reached by modeling lunar occultation light curves for molecular 
contamination effects. 

The two SR variables rj Gem and Gem clearly do not show any variation in their angular 
diameter in K and L/ bands unlike the Miras. SW Vir and ji Gem appear to be candidates for 
overtone mode from their Q-values. However as we are using occultation UD radii which are greater 
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than Rosseland radii, our Q-values may be less than the theoretical values. The case of r] Gem is 
complicated by the presence of two periods in its optical light curve. Further high angular resolution 
studies of SR variables are clearly warranted. 

The authors wish to thank the referee for many valuable suggestions for improving the paper. 
This research has made use of the AFOEV and SIMBAD databases, operated at CDS, Strasbourg, 
France. This work was supported by Dept. of Space, Govt, of India. 
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Fig. 1. — Occultation light curves of U Ari in the K-band: model fit, residuals of fit (lower panel) 
and convergence of fit (inset). The dotted and solid line is the observed data and model fit curve. 
The bestfit uniform disk (UD) angular diameter is 7.3±0.3 mas. 
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Fig. 2. — Occupation light curves of Z Sco in the K-band: model fit, residuals of fit (lower panel) 
and convergence of fit (inset). The dotted and solid line is the observed data and model fit curve. 
The bestfit uniform disk (UD) angular diameter is 3.8±1.0 mas. 
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Fig. 3. — Occultation light curves of r\ Gem on 8 January 2001 (top) and 4 March 2001 (bottom) 
and : model fit, residuals of fit (lower panels) and convergence of fit (insets) in the K-band (left) 
and L'-band (right). The dotted and solid line is the observed data and model fit curve. The 
best-fit uniform disk diameters are 12.7±0.3 mas at K and 12.7±1.0 mas at V (8 Jan 2001) and 
12.8±0.3 mas at K and 12.8±2.0 mas at L' (4 Mar 2001) respectively. 
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Fig. 4. — Occultation light curves of fx Gem in the K-band (left) and L' -band (right) : model fit, 
residuals of fit (lower panels) and convergence of fit (insets). The dotted and solid line is the 
observed data and model fit curve. The best-fit uniform disk angular diameter is 13.7±0.5 mas at 
K-band and 14-8±1.0 mas at L '-band. 
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Fig. 5. — Occultation light curves of SW Vir in the K-band: model fit, residuals of fit (lower panel) 
and convergence of fit (inset). The dotted and solid line is the observed data and model fit curve. 
The bestfit UD angular diameter is 15.9±0.6 mas 
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Fig. 6. — Occultation light curves of TV Gem in the K-band : model fit and convergence of fit 
(inset). The dotted and solid line is the observed data and model fit curve. The best-fit UD angular 
diameter is 4-8±0.2 mas. Lower panel of the graph shows the enlarged residuals particularly on 1st 
and 2nd fringe positions which is marked by arrow for star+shell model fit (solid line) and single 
star model fit (dotted line). 
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Fig. 7. — The blackbody fit to photometry data available in the literature for U Ari, U Ori, Z Sco 
and r\ Gem. 
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Fig. 8. — Linear radius vs. period plot for the Miras and semi-regular variables . The filled circles 
represent our LO measurements in the K-band. The fundamental and first- overtone curves in the 
mass range 1.0 Mq -2.0 Mq are derived from the radius-period relation of eqn. 3 and eqn. 4- 
Some of available linear radii for U Ori and U Ari in the literuare are superimposed. 
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Table 1: Log of Observations 



Star (var.) 


IRC 


Date of 


Filter(/xm) 


PA 




event type a 






obser. 


(A/AA) 


(deg.) 


(km s _1 ) 




U Ari (Mira) 


+10040 


19 Feb 2002 


2.2/0.4 


18 


0.475 


D 








3.8/0.6 








Z Sco (Mira) 


-20306 


22 Mar 2003 


2.2/0.4 


155 


0.651 


R 


U Ori (Mira) 


+20127 


13 Mar 2000 


2.2/0.4 


136 


0.66 


D 


SW Vir (SRb) 


+00230 


01 Jun 2001 


2.2/0.4 


74 


0.504 


D 


r] Gem (SRa) 


+20139 


08 Jan 2001 


2.2/0.4 


21 


0.330 


D 








3.8/0.6 












04 Mar 2001 


2.2/0.4 


68 


0.645 


D 








3.8/0.6 








jj, Gem (Lb) 


+20144 


04 Mar 2001 


2.2/0.4 


77 


0.767 


D 








3.8/0.6 








TV Gem (SRc) 


+20134 


14 Nov 2000 


2.2/0.4 


212 


0.437 


R 



a D: Disappearance; R: Reappearance 
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Table 2: Individual source parameters 



Parameters U Ari Z Sco U Ori rj Gem /j, Gem SW Vir TV Gem 

Spectral type 1 M4-9.5 Hie M4/5 Hie M6.5-9 Hie M3 III M3 III M7 III Ml-0 lab 

V mag 1 7.2-15.2 8.7-13.4 5.3-12.6 3.20-3.90 3.20 8.2-9.4 8.7-9.5 

(max-min) 

K mag 2 1.59 1.48 -0.49 -1.49 -1.89 -1.87 0.99 

L' mag 2 0.28 1.50 -1.41 -1.59 -2.01 -2.28 0.62 

Variability type Mira Mira Mira SRa Lb SRb SRc 

Period (days) 1 371 352 372 233 27 10 150 182 

Distance (pc) 776±155 3 778±145 3 306±61 3 107±22 4 71±5 4 143±24 4 1200±300 5 

Luminosity 11 8.4±3.5 8.3 ±3.3 9.7±4.0 2.3±1.0 1.8±0.36 4.7±1.7 68.5±34.9 
(1(P L Q ) 

Mass loss 57 6 - 29 7 1.4 8 0.44 s 17 7 200 9 
(10~ 8 M Q yr- 1 ) 

Outflow velocity 6.0 6 - 7.5 7 - - 7.8 7 12 9 

(km s^ 1 ) 

Table References: 1. SIMBAD database; 2. Gezari et al. 1999; 3. Whitelock et al. 2000b; 4. Pcrryman 

et al. 1997 (Hipparcos catalog); 5. Underbill 1984; 6. Winters et al. 2003; 7. Knapp et al. 1998; 8. Drake 

et al. 1986; 9. Loup et al. 1993; 10. Percy et al. 2001. 



li 



Luminosity is estimated using the relation L = 47rd* 2 Fb ; . 
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Table 3: Derived Infrared angular diameters 



btar 


Date 


Phase 


\ / A \ 

A/ AA 

w 


Ang. JJia (UDJ 
(mas) 


U Ari 


19 Feb 2002 


0.57 


2.2/0.4 


7.3±0.3 


Z Sco 


22 Mar 2003 


0.26 


2.2/0.4 


3.8±1.0 


U Ori 


13 Mar 2000 


0.33 


2.2/0.4 


11.9±0.3 


SW Vir 


01 Jun 2001 


0.56 


2.2/0.4 


15.9±0.6 


r] Gem 


08 Jan 2001 




2.2/0.4 


12.7±0.3 




08 Jan 2001 




3.8/0.6 


12.7±1.0 




04 Mar 2001 




2.2/0.4 


12.8±0.3 




04 Mar 2001 




3.8/0.6 


12.8±2.0 


\i Gem 


04 Mar 2001 




2.2/0.4 


13.7±0.5 




04 Mar 2001 




3.8/0.6 


14.8±1.0 


TV Gem 


14 Nov 2000 




2.2/0.4 


4.8±0.2 
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Table 4: Previous Angular size measurements from Optical to Near-IR 



Star 


Date of Obs. 


Method" 


Phase 


A/AA 


Ang. Dia(UD) 


Reference 










(/an) 


(mas) 




U Ari 


03 Sep 1977 


LO 


0.49 


1.62/0.42 


6.11 ±0.34 


Ridgway et al. 1979 


(Mira) 














U Ori 


15 Oct 2000 


LBI 


0.88 


2.20/0.40 


15.59±0.06 


Mennesson et al. 2002 


(Mira) 


15 Nov 2000 


LBI 


0.96 


3.75/0.70 


25.66±0.69 


Mennesson et al. 2002 




Oct 2000 


LBI 


0.83 


2.20/0.1 


10.6 


Perrin et al. 2004 




Nov 2001 


LBI 


0.91 


2.20/0.1 


9.66±0.12 


Perrin et al. 2004 


SW Vir 


26 Jan 1981 


LO 


0.02 


1.62/0.04 


16.82 ±0.34 


Ridgway et al. 1982 


(SRb) 


01 Sep 1981 


LO 


0.48 


2.17/0.03 


16.11±0.13 


Schmidtke et al. 1986 




29 Jun 1982 


LO 


0.48 


2.28/0.40 


16.77±0.23 


Schmidtke et al. 1986 




29 Feb 2000 


LBI 


0.51 


2.20/0.40 


16.24±0.06 


Mennesson et al. 2002 




12 Mar 2000 


LBI 


0.56 


3.75/0.70 


22.88±0.33 


Mennesson et al. 2002 


rj Gem 


1988-90 


LBI 




0.55/0.02 


11.43±0.55 


Mozurkewich et al. 2003 


(SRa) 


1988-90 


LBI 




0.80/0.02 


10.91±0.11 


Mozurkewich et al. 2003 




1993 


LBI 




0.712/0.012 


11.75±0.27 


Quirrenbach et al. 1993 




1993 


LBI 




0.754/0.005 


10.70±0.15 


Quirrenbach et al. 1993 




31 Mar 2001 


LO 




2.20/0.40 


12.57±0.04 


Richichi et al. 2003 


\x Gem 


1988-90 


LBI 




0.55/0.02 


13.48±0.19 


Mozurkewich et al. 2003 


(Lb) 


1988-90 


LBI 




0.80/0.02 


13.99±0.14 


Mozurkewich et al. 2003 




1993 


LBI 




0.754/0.005 


13.50±0.13 


Quirrenbach et al. 1993 




1993 


LBI 




0.712/0.012 


13.97±0.28 


Quirrenbach et al. 1993 




1987 


LBI 




2.2/0.40 


13.50±0.15 


Di Benedetto et al. 1987 


TV Gem 


15 Aug 1982 


LO 




0.55 


5.31±0.91 


Radick et al. 1984 


(SRc) 


30 Mar 1993 


LO 




2.2/0.40 


4.9±0.30 


Ragland et al. 1997 




03 Feb 1993 


LO 




2.2/0.40 


4.46±0.07 


Richichi et al. 1998 



a LO: Lunar Occultation; LBI: Long Baseline Interferometry 
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Table 5: Derived angular diameters, bolometric fluxes and effective temperatures 



Star 


Spt.Type 


mfc 


m„ 


UD (K-band) 
mas 


F b oi xlO 8 
(erg cm -2 s _1 ) 


T e// 
(K) 


U Ari 


M9 Hie 


1.26 


14.50 


7.3±0.3 


45±5 


2250±80 


Z Sco 


M4/5IIIe 


1.33 


11.6 


3.8±1.0 


46±5 


3120±420 


U Ori 


M8 Hie 


-0.88 


10.89 


11.9±0.3 a 


336±33 


2905±80 


SW Vir 


M7 III 


-1.74 


7.90 


15.9±0.6 


735±110 


3060±130 


i] Gem 


M2.5 III 


-1.49 


3.70 


12.8±0.3 


760±105 


3450±125 


fi Gem 


M3 III 


-1.89 


3.20 


13.7±0.5 


1140±170 


3675±140 


TV Gem 


Ml la 


1.16 


6.83 


4.8±0.2 


153±15 


3750±120 



a The angular diameter of U Ori is taken from Mondal & Chandrasekhar (2004) 



-31 - 



Table 6: Linear radii and distances of Mira and SRs 



Star 


K-band UD 

(mas) 


Period 

(days) 


PL dist. 

(PC) 


Hip. dist. 

(PC) 


adopted dist. 

(PC) 


Lin. radii 

(R©) 


Q- value 
(for 1 M ) 


U Ari 


7.3±0.3 


371 


776±155 a 




776±155 


610±125 


0.025 


Z Sco 


3.8±1.0 


352 


778±145 a 


585±875 b 


778±145 


320±100 


0.062 


U Ori 


11.9±0.3 


372 


306±61 a 


505±420 b 


306±61 


391±78 


0.048 


SW Vir 


15.9±0.6 


150 


98±19 a 


143±24 


143±24 


244±42 


0.039 


r\ Gem 


12.8±0.3 


233 
20 c 


150±30 a 


107±22 


107±22 


146±30 


0.130 
0.012 


fi Gem 


13.7±0.5 


27 




71±5 


71±5 


104±8 


0.025 


TV Gem 


4.8±0.2 


182 




1492±2340 


1200±300 d 


623±158 


0.012 



a From Period- Luminosity (PL) relation of Whitelock & Feast (2000b). 
b Revised Hippaecos parallax from Knapp et al. (2003). 

c This secondary period is determined from visual photometric observations by Percy et al. (2001). 
d Adopted from Underhill (1984). 



